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The Hamiltonian for a system of relativistic bodies interacting by their gravitational 
field is found in the post-Minkowskian approximation, including all terms linear in the 
gravitational constant. It is given in a surprisingly simple closed form as a function of 
' canonical variables describing the bodies only. The field is eliminated by solving inhomo- 

geneous wave equations, applying transverse-traceless projections, and using the Routh 
functional. By including all special relativistic effects our Hamiltonian extends the results 
described in classical textbooks of theoretical physics. As an application, the scattering 
of relativistic objects is considered. 

^ '. 

{^j^' Keywords: Post-Minkowskian approximation 

Post-Newtonian (PN) approximation methods in general relativity are based on 
the weak-field limit in which the metric is close to the Minkowski metric and the 
assumption that the typical velocity t; in a system divided by the speed of light is 
very small. In post-Minkowskian (PM) approximation methods only the weakness 
of the gravitational field is assumed but no assumption about slowness of motion 



' is made. In the PM approximation we obtain the Hamiltonian for gravitationally 



interacting particles that includes all terms linear in gravitational constant G. It 
' thus yields PN approximations to any order in 1/c when terms linear in G are 

. considered; and it can also describe particles with ultrarelativistic velocities or with 

zero rest mass. 

We use the canonical formalism of Arnowitt, Deser, and Misner (ADM)2J where 
. the independent degrees of freedom of the gravitational field are described by hf^^, 

?H ' the transverse-traceless part of hij = gij —Sij [h^^^ — 0, hj^, = 0, z, j — 1, 2, 3), and 



by conjugate momenta / {16ttG)tt^^ '^'^ . The field is generated by N particles with 
rest masses ma located at x^, a — 1, ...N, and with momenta Pa- We start with the 
HamiltoniarP correct up to G^ found by the expansion of the Einstein equations 
(the energy and momentum constraints) in powers of G and by the use of suitable 
regularization procedures. When we consider only terms linear in G and put c — 1 
this Hamiltonian reads 

i/n„=E^-|GE"fl + 4 + 4) (1) 
^ 2 Tab V mi mf^J 



7G V — (7pa-p6 + (Pa-na6)(pb-nab)) - ^ V h. 

^ Tab 2 ^ rrir 
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where rua — (m^ + p^) ^ , Uabrab — Xa — Xb, |nafc| — 1. The equations of motion 
for particles are standard Hamihon equations, the Hamihon equations for the field 
read 

■ijTT Tr r< xTTij -rpj, iR^^ cTTkl 

TT^ = -167rG(5;,; ^ = 167rG<5y J^^mtt ^ (2) 

here the variational derivatives and the TT-projection operator (5^;"^*"' = 
i {AikAji + AuAjk - Ay Afe;) A"2, Ay = %A - di dj, appear. These equations 
imply the equations for the gravitational field in the first PM approximation to 
be the wave equations with point-like sources ^ S'-^^x — Xq). Since both the field 
and the accelerations Pa are proportional to G, the changes of the field due to the 
accelerations of particles are of the order 0{G^). Thus, in this approximation, wave 
equations can be solved assuming field to be generated by unaccelerated motion of 
particles, i.e., it can be written as a sum of boosted static spherical fields: 

/^rW = E ^ 1 ' I ; ^"^^ . (3) 

b ' ' W f — x^ sm ^6 

where x — Xq = riajx — Xaj and cos^a = ^a^a/\^a\- Surprisingly, it is possible to 
convert the projection Sfj^ (which involves solving two Poisson equations) into an 
inhomogeneous linear second order ordinary differential equation and write 

/^rW = E n! [yP' - (n''-P'')'(32/ + 2)] S.J (4) 

\x - xb\ 77ib y(l + y)- ^ 

+ 2 [1 - x2(l - 2cos2 9b)] pMPbj + [(2 + y) {nb-pb)'^-{2 + 3y - 2x1) pI] nb.nbj 

+ 2(n6-pb) + 2y) {nuPbo + Pbtnbj) } + 0{mb±b - Pb)G+0{G^) ; 



here y = yb — y 1 — x^ sin^ 0h and we anticipate 0{mbXb — Pb) G. 
In the next step we use the Routh functional (see, e.g., Ref.!^ 

i?(x„. Pa, hl^, hl^) ^H^ ^ j CPX TT^^^^ hj^ , (5) 

which is "the Hamiltonian for the particles but the Lagrangian for the field." Since 
the functional derivatives of Routhian vanish if the field equations hold, the (non- 
radiative) solution ([3]) can be substituted into the Routh functional without chang- 
ing the Hamilton equations for the particles. Using the Gauss's law, an integration 
by parts and similar standard steps (such as dropping out total time derivatives, 
i.e. a canonical transformation) and the explicit substitution for hjj" [x = x^) we 
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get the Hamiltonian for a A'^-particle gravitating system in the PM approximation: 

Hlin = V'^a + tG — (7pa-p6 + (Pa-nafc) (pfe -nab)) - Ilhl—l 

4 Tab 2 ^ Tab 

1 (maTOfc)"^ 



X 1 + 



pI 



pI 



4 ^ Tab (j/ba + l)^2/ba 
a,b^a 



2 2(pa-Pb) (Pb-nba) 



(6) 



-2(pa-nbQ)(pb-nba)(Pa-Pb)Pb + (Pa-nba) Pb - (PaPb) Pb =2 + ^ -Pa(Pb-nba) 



+ (pQ-nba)^(pb-nbQ)^ + 2(pa-nba)(pb-nba)(Pa-Pb) + (Pa-Pb)^ - (pQ-nba)^Pb 
- 3p^(pb-nbQ)^ + (pa-nba)^(Pb-nba)^ + 8(pa-nba)(Pb-nba)(Pa-Pb) 



PaPb - 3(pa-nba)^P^ 



Uba 



Vba 



nib 



— V ml + (nba-Pb) 



Since the PM approximation can describe uhrarelativistic or zero-rest-mass par- 
ticles, we calculated gravitational scattering of two such particles using the Hamilto- 
nian ([6]). If perpendicular separation b of trajectories (|b| is the impact parameter) 
in the center-of-mass system (pi — — p2 = p) is used, p-b — 0, we find, after evalu- 
ating a few simple integrals, that the exchanged momentum in the system is given 

by 



Ap 



b^ IpI mi + m2 



1 



1 



+ 



TO 17772 



(7) 



The quartic term is all that remains from the field part hjj in agreement with 
WestpfahP who used a very different approach. 

The Hamiltonian ([6]) can also describe a binary system with one massless and 
one massive particle orbiting around each other. This is not obvious: the second, 
fourth or even sixth-order PM approximation would not be able to describe massless 
test particles orbiting around a Schwarzschild black hole. 
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